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The specification has been amended and a substitute 
specification is attached (both a clean and marked up 
copy) . The drawings have been amended and three sheets of 
replacement drawings are attached. 

In the claims, claims 2, 4, 11 and 12 have been 
amended. For example, the term "basic" has been changed to 
"base." The term antialias has been changed to 
"antialiasing." The term "amount" has been changed to 
"magnitude" and the term "cycle reduced" has been changed 
to "downsampled . " 



No new matter has been added. 
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SUBSTITUTE SPECIFICATION 
METHOD FOR THE LINEARIZATION OF FMCW RADAR DEVICES 

CROSS REFERENCE TO RELATED APPLICATIONS 
This application is based upon German Application 
Serial No. 10315012.9 Filed on April 2, 2003 wherein 
priority is claimed under 35 U.S.C. 119. 

BACKGROUND OF THE INVENTION 

The invention relates to a method for the linearization 
of a frequency modulated continuous wave (FMCW) radar 
device . 

With FMCW radar devices, phase errors occur in the 
reception signal, as the result of a not strictly linear 
incrcaoc characteristic in a transmission frequency of an 
FMCW radar device. This results in a loss of accuracy in 
the evaluation of the signal. 

In this case, various methods are known for linearizing 
an FMCW radar device. For example, one known method 
discloses a method wherein a delayed transmission signal is 
generated via an optical delay line having a defined running 
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time. In this case, the delayed transmission signal is 
mixed with the non-delayed transmission signal. Thus, it is 
possible to determine the deviation from the ideal linear 
characteristics. In particular, the disadvantages of this 
method are the great amount of construction effort, expense 
and the age-related clouding of the optical components. 

In addition, it is also known to generate a signal 
function that can be predetermined, such as a linear 
characteristic via a digital direct synthesis (DDS) system. 
However, a disadvantage of this system is that there is a 
large construction effort and expense associated with this 
cost intensive structure. 

An FMCW radar system having linear frequency modulation 
is known from German Patent DE 4104792 Al . In this system, 
there is an oscillator that generates a ramp- shaped 
modulated signal having a modulated oscillation frequency, 
which is emitted by an antenna. A monitoring and control 
unit determines the deviation of the measured modulated 
oscillator frequency of a signal reflected from an object 



2 



via a Hilbert transformation using a linearized reference 
value. In this case, a phase deviation can be determined 
from the deviation that has been determined. 

From this deviation, a correction signal is generated, 
which is passed to the oscillator, whereby the oscillation 
frequency is modified in accordance with the predetermined 
reference values. A disadvantage of this system is that 
this involves a great construction effort, which is cost 
intensive, since additional components are required. 

SUMMARY OF THE INVENTION 

The invention relates to a method with which an FMCW 
radar device can be linearized, without a great construction 
effort or expense. 

In this case, a correction phase term is calculated on 
the receiver side for compensation of the phase error in the 
reception signal q(t) . In contrast with known methods, no 
costly components are required in the transmission branch of 
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the FMCW radar device when using this method and this thus 
results in reducing costs. 

This method, which relates to the invention should 
include the following steps: 

selecting a number L of consecutive ramp- shaped 
reception sequences q k (n) of the reception signal, which 
number can be predetermined, with K=l, . . . ,L, the phases 
arg{q k (n)} of which can be represented as a polynomial of the 
N th order for the time index n, with the polynomial 
coefficient m c , with 0=1,..., N, 

transforming a spectrum range Q(e jn )of the selected 
reception sequences q(n) into a basic band reception 
sequence q k (ji) with k=0,...L-l, 

iteratively calculating a correction phase term for 
partial compensation of non-linear frequency components in 
the basic band reception sequences q k {n) by means of 
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calculating polynomial coefficients of the individual 
basic band reception sequences q k (n) via estimating methods 

wherein ^ki n ) are the sequences that have already been 
iteratively phase-corrected, and whereby the iteration is 
stopped once the parameter change between two consecutive 
iterations, which can be predetermined, remains below a 
threshold e, which can be predetermined. 

In another embodiment of the invention, the scanning 
cycle T A of the Kffl ramp signals is reduced, in a method step 
after the baoic base band transformation. For this purpose, 
a prior Antiliao antialiasing low pass filtering of the ramp 
signals q k (n) is performed. In addition, it is beneficial if 
the downsampling factor K, by which the scanning cycle T A is 
reduced, lies between K=30 and K= 60. 

In contrast, with an alternative linearization of the 
frequency modulation of the radar transmission signal, the 
linearization is performed on the receiver side, via a 
digital processing of the reception signal. Thus, no 
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additional components are required for this method. The 
calculations required for this method can be performed by 
the components for signal evaluation that are already 
present. This results in additional advantages with regard 
to the production costs. 

Another advantage of the invention is that there can be 
in an inexpensive estimation method that can be easily 
implemented, since the reception signal is processed in a 
the baoic base band. This signal results in a reduction of 
the scanning cycle T A by a factor K between 30 and 60. In 
addition another advantage of this method is that there can 
be a calculation of the polynomial coefficients via a robust 
estimation method. This is even in the case of channels 
subject to severe interference. Thus, there can be a 
relatively good compensation of phase errors with signal to 
noise ratios (SNR) of up to -20 dB using this method. Thus, 
the method can be used for not only the civil sector, but 
also in the military sector or for distance measurements. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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Other objects and features of the present invention 
will become apparent from the following detailed description 
considered in connection with the accompanying drawings, 
which discloses at least one embodiment of the present 
invention. It should be understood, however, that the 
drawings are designed for the purpose of illustration only 
and not as a definition of the limits of the invention. 

In the drawings, wherein similar reference characters 
denote similar elements throughout the several views: 

FIG. 1 is an exemplary idealized ramp- shaped FMCW 
transmission signal; 

FIG. 2 is an exemplary representation of the refection 
of the radar signal on a point shaped object at distance R; 

FIG. 3 is a schematic diagram of the method that 
relates to the invention; 
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FIG. 4 is a representation of a baoic the base band 
shift of a selected spectrum range in the reception signal; 

FIG. 5a is a representation of a spectrum of a cycle 
reduced downsampled basic band signal before and after the 
band pass filtering; and 

FIG. 5b is a representation of an amount a magnitude 
spectrum of a radar reception signal before and after the 
phase correction. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Turning now in detail to the drawings, FIG. 1 shows the 
starting point of this method is the transmission signal of 
an FMCW radar that is in the general polynomial approach of 
the N th order for a transmission ramp as: 

r ^ ( r mi 2 mi 3 mi 4 V 
x(t) = zxv{j27\j 0 t + —t 2 + — t l +— t 4 +...+y 0 J} (1) 

With t e[0,T r ] . With this design, the coefficients m c 
with 0=2, . . . ,N represent the polynomial parameters, y° 
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represents the initial phase, and T r represents the time 
duration of the transmission ramp. The parameter mi or slope 
that is assumed to be given, indicates the linear component 
mi=B/t, which is exclusively present in an ideal case, 
wherein B describes the frequency otrokc deviation of the 
FMCW radar system. Thus in the following formula, a single 
reflecting point shaped object is assumed. Thus, the 
corresponding parameters are represented in FIG. 1, using 
the example of an idealized transmission signal. Thus with 
Equation (1) the reception signal in an ideal, noise free 
case is then: 



Since the determination of relative distances between 
objects that are moving relative to the radar sensor at 
approximately the same speed represents a possible use of 
the method described, it can be assumed that s(t) = 2R in 
the following (Fig. 2) . 



A(l). 



s(t) 



} (2) 
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For the momentary wavelength X(t), the following 
results at the speed of light c: 



^-_m=iU + ± m A wi,h m --^m--n+± m / 

Mf) c cV tl ) In dt M 

It is advantageous if the reception signal y(t) is mixed 
into a lower frequency position. The mixed reception signal 
q e (t) results as follows, from Equations (1) and (2): 



q e (t) = y{t) x (0 = exp - j—\ fo + L m t t 

If q e (t) describes a harmonic exponential function, i.e. 
if m c ^ 0 with { = 2, . . . , N, then the ideal case is present, 
and the transmission signal x(t) possesses a purely linear 
increase in the transmission frequency. 

In order to indicate the maximal deviation from this 
ideal transmission frequency progression, a maximal relative 
linear error is defined as 
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\m-f ideal (t)\ 

0:= max- — 



i 

= — max 
/ 0 -0 m, '4 0 ^1 



*=2 



with fid(t) =f 0 +mit for the ideal case. The symbol 6 
therefore describes a maximum deviation of the momentary 
frequency progression for f 0 = 0 from the linear function 

f id {t) q . Taking into consideration the noise term r(t) that 

is always present in real applications, the following is 
obtained from Equation 3, with q(t) = q e (t)+r(t) 



qit> exp^ - j 



/o + Z'M'H + KO (4) 



If the complex-value noise signal r(t) is assumed to be 
white and gauoo gaussian distributed in Equation 4, for the 
sake of simplicity, the noise influence for a significantly 
large signal-to-noise diotancc ratio can also be considered 

as phase noise e jr{t) , whereby r(t) is real -valued, white, and 

gauoo gaussian distributed. The phase signal, which is 
made time-discrete with the scanning cycle T A , results from 
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arg 



{q{nT A )\ = 4(nT A 



)=-^(/ 0 + Em € (nrjM+;(nr il ) 



(5) 



Proceeding from the known (received) phase scanning 
values <|>(nT A ), the method according to the invention 
implements a robust estimate of the polynomial parameters m P/ 
with ? = 2, . . . , N. In an advantageous embodiment of the 
invention, a robust estimate of the unknown object distance 
R is additionally possible. 

Fig. 3 shows a block schematic in the form of a 
sequence schematic of the method according to the invention. 
The initialization 1 takes place, according to the 
invention, by means of a baoic base band shift BV of L 
consecutive ramps of the sequence q(n) = q(nT A ) that results 

from Equation (4) , by means of scanning time discretitation . 
The shift of a spectrum range, which can be predetermined, 

of Q( e ^) = F {q( n )} into the baoic base band is shown as an 

example in Fig. 4. The cut-off frequencies Q x and Q u are 
defined, for example, by way of the minimum and maximum 
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distances R min and R max , which are determined by the purpose 
of use of the FMCW radar system: 

= —R n ^T A , Q„ = —R^m,T A , Q,„ = l(fl. - Q,) 
c c i 

The ramp signals q k (n) that are formed in this manner 
possess a band-pass spectrum. It is advantageous that after 
the baaic base band shift and the subsequent Antialiao 
antialiasing low-pass filtering TF, the scanning cycle can 
be reduced by a factor K. In this connection, the factor K 
particularly lies between K = 30 and K = 60. This results 
in advantages with regard to a reduction in effort/expense 
of all further partial operations. 

The next step or progression of the block diagram in 
Fig. 3 shows the sequence of the iterations for the 
calculation of a correction phase term for partial 
compensation of non- linear frequency components in the basic 
band reception sequences qf{m) with m as a time index after 

cycle reduction. In the further process steps, a signal 
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sequence whose scanning cycle T A was reduced is assumed with 
qf{m) . Of course, these process steps can also be carried 
out with a sequence qf{n) f whose scanning cycle T A was not 
reduced. 

Within the loop that runs through in each iteration, it 
is advantageous if band-pass filtering of the sequences 

qf{m) is located at the beginning, which is advantageously 
carried out as window formation windowing in the spectrum 
range, whereby all the spectral coefficients outside the 
window are set to zero. It is advantageous if a Hamming 
window is used as the window function. 

In the case of several reflecting objects, the 
reception signal spectrum of a single object is selected by 
means of the band-pass filtering. All of the other, 
undesired spectra are suppressed. In addition, suppression 
of those noise-like spectrum components that do not belong 
to the working component takes place. In this way, the SNR 
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can be increased after the band-pass filtering, even in the 
case of a single reflecting object. 

To calculate the band-pass filtering, it is 
advantageous to perform a faat Fast Fourier tranof ormation 
Transformation FFT of the sequences q[ n (m) - For this 
purpose, it is advantageous if the individual discrete 
Fourier tranoformcd transforms Q^iju) of the basic base band 

reception sequences qf(m) are calculated, whereby 

$\p) = FFT[qf(m)} for k = 1,..., L. The point count M of the 
FF¥ number of FFT points M typically lies at up to 1000. In 
the case of a reduction of the scanning cycle T A that has 
advantageously been performed, and prior Antialiao 
antialiasing low-pass filtering, the FFT can be carried out 
with a longer point count smaller number of points M of 64 
to 256. 

In a next iteration step, filtered baaic base band 
reception sequences q { k l) {m) are calculated by means of a band- 
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pass filter according to Q k {i \ju) = Hm) Ql° (M) < with 
(% s \ M ) = FFT{qj; 0 (m)} ' whereb Y is a spectrum window that can 

be predetermined, and fie[fi m9 fi t ] indicates the range of the 
spectrum window having a window center point |a max that can be 
predetermined, whereby, with a lower limit \i u that can be 

predetermined, and an upper limit |i i that can be 
predetermined. In order to determine the position of the 
center point \i max/ the FFT amount magnitude spectra Q^\/j) 

are averaged over the number L, and the maximum is sought 
from the averaged spectrum, according to the following 
formula : 



G (0 (/O = max -Vlel 0 ^) 



For example, for a single reflecting object, \i 1 = 0 and 
(i u = M-l. For several non-overlapping objects, the limits jj. 1 
and n u must be selected so that in each instance, only the 
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spectral component that corresponds to the object of 
interest is detected. 

Subsequently, it is advantageous if the calculation of 
the individual inverse discrete Fourier transformed q^{m) of 

the filtered FFT spectra Q^(ju) takes place, whereby 
q { k n (m) = IFFT{Q^(M)} for k = X ' ■ • L ' 

In a subsequent iteration step, the distances Rk (l) are 
estimated by means of a maximum likelihood estimating 
method. The polynomial coefficients Rk {l) can be estimated 

from the estimated distances ™e,l . Subsequently, the 

distances R k are averaged over L reception sequences q k {ri) . 

In a final iteration step, a reception sequence q { * l) (m) with 

the averaged estimated polynomial coefficients ™/ is 
calculated as the starting point for a next iteration loop, 
according to 
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The maximum likelihood parameter estimate takes place by 
using the linear estimator 

v = [A T A\~ l A T b (6) 

with the vectors v = [Rf Q9 Rm l9 Rm 29 ...] T , b = [®(0)MT A ), <&{2T A ) 9 ..] T and 

a matrix A having the dimension L q x (N+l) , whereby L q 
represents the length of a sequence q^{m) - The 

determination of the phase values <t>(nT A ) in the vector b 
takes place by means of a phase unwrapping method. 

A lcaocr reduced linearity error in the cycle reduced 
beHEHre downsampled base band signal sequence q { k l) (m) is 

successively achieved by means of the phase correction, and 
this results in a lesser estimate variance, in each 

instance. If the condition - R u) \ < € has been fulfilled 
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after i = I iterations, final estimated values R and m e 

can be calculated from the individual intermediate results 

R (i) and , which have been averaged over all L reception 
sequences in the i th iteration, in each instance, according 
to : 

R = R {,) and m ( -it^j 0 , *-2,...N 

& 1 = 1 

In the following, it will be shown that it is possible, 
using the method according to the invention, to reconstruct 
the spectral resolution, which has been reduced as a result 
of the linearity error in the FMCW radar transmission 
signal, of the ocannod sampled , mixed reception signal q(n) . 
This will be shown in exemplary manner below, using the 
example of two objects at different distances. 

The starting point is a quadratic transmission signal 
sequence (chirp) with a linearity error of 6=5%, whereby 
the simulation parameters are selected to be N = 2, m 2 = 
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5xl0 12 , SNR = -18 dB, K = 40, and M = 64 . For the other 
system parameters, the following table applies: 



Parameter 


Value 


Minimum possible distance 


Rmin = 900 m 


Maximum possible distance 


Rmax = 1,100 m 


Ramp increase time 


T r = 1 ms 


Daoic Radio Fundamental 
frequency 


f 0 = 35 GHz 


Frequency stroke deviation 


B = 100 MHz 


Scanning Sampling frequency 


f A = 2 MHz 


Spectrum averaging L 


L = 20 


Spectrum averaging, total 


25 x L = 500 


Spectrum window w ( \i) 


Hamming window 



The objects are located at the distances R 0 = 900 m and Ri = 



1,100 m. Fig. 5 shows the amount magnitude spectra 

of the cycle reduced downsampled signal q[ 0) (m) 

for a fixed k = 0, . .., L-l during the first iteration 
before (top representation in Fig. 5a) and after (bottom 
representation in Fig. 5b) band-pass filtering. The partial 
spectra, which can be easily separated, are clearly evident, 
whereby the bottom partial spectrum, which belongs to the 
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object at the distance R lt is selected in the band-pass 
filtering . 



The top figure in Fig. 5b shows the range of the amount 
magnitude spectrum of the radar reception signal that is 

of interest, which is established by the minimum and maximum 
distance (see table) . As the result of the non-linearity in 
the transmission signal, the frequency lines that correspond 
to the two objects are smeared and can hardly be recognized. 
Using the parameters estimated for the object at the 
distance R lt using the method according to the invention, the 
inherently non-linear component of this signal is now 
compensated, as shown in the bottom part of Fig. 5b. 

A comparison of the figures shows the clear increase in the 
object resolution. The small triangles in the bottom part 
of Fig. 5b characterize the positions of the frequency lines 

in the amount mangitude spectrum \Q{ eJn ) , which would generate 
a transmission signal having a linear frequency increase. 
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Accordingly, while at least one embodiment of the 
present invention have been shown and described, it is to be 
understood that many changes and modifications may be made 
thereunto without departing from the spirit and scope of the 
invention as defined in the appended claims. 
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ABSTRACT OF THE DISCLOSURE 

A method for the linearization of frequency modulated 
continuous wave (FMCW) radar devices having non- linear, 
ramp-shaped, modulated transmitter frequency progression 
x(t) . With this invention, a correction phase term for 
compensation of the phase error in the reception signal q(t) 
is calculated on the receiver side in this device. 
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